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Electromagnetic interactions between colliding heavy ions at the Large Hadron Collider (LHC) 
at CERN will give rise to localized beam losses that may quench superconducting magnets, apart 
from contributing significantly to the luminosity decay. To quantify their impact on the operation 
of the collider, we have used a three-step simulation approach, which consists of optical tracking, 
a Monte-Carlo shower simulation and a thermal network model of the heat flow inside a magnet. 
We present simulation results for the case of 2t)8 Pb 82 + ion operation in the LHC, with focus on the 
ALICE interaction region, and show that the expected heat load during nominal 208 pb 82+ operation 
is 40% above the quench level. This limits the maximum achievable luminosity. Furthermore, we 
discuss methods of monitoring the losses and possible ways to alleviate their effect. 

PACS numbers: 29.20.db, 25.75.-q, 84.71.Ba 



I. INTRODUCTION 

The Large Hadron Collider (LHC) is presently being 
commissioned at CERN pQ. Its main beam parameters 
are given in Table [TJ In its first phase of operation it will 
collide protons but later also heavy nuclei, starting with 
208 Pb 82+ at center-of-mass energies up to 1.15 PeV [2]. 
This will open up a new energy regime in experimen- 
tal nuclear physics, extending the study of the hadronic 
matter ( "quark-gluon plasma") that existed in the early 
universe about lCP 6 s after the Big Bang [3J. At the 
same time, colliding heavy-ion beams at this unprece- 
dented energy will present beam physics challenges not 
encountered in previous colliders. 

When two fully stripped ions collide at an interaction 
point (IP), a variety of processes leading to fragmenta- 
tion and particle production can occur. Hadronic nuclear 
interactions, which are usually the main object of study 
of the experiments, occur only when the impact param- 
eter b is smaller than about twice the nuclear radius R. 
Ultraperipheral collisions are those in which two collid- 
ing ions pass close to each other with b > 2R. In such 
events, the intense Lorentz-contracted fields of the nu- 
clei can be represented as pulses of virtual photons in 
the equivalent photon picture of Fermi, Weizsacker and 
Williams [4]. These extremely energetic photons collide 
and cause electromagnetic interactions that are particu- 
larly strong in heavy-ion collisions because the density of 
virtual photons around the nuclei is proportional to Z 2 . 
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TABLE I: Design parameters for the LHC's proton and 
208 Pb 82+ beams in collision conditions [J. The values of C 
and (3* refer to IP2 for 208 Pb 82+ ions and IP1 and IP5 for 
protons (see Fig. [TJ . 



Particle 


P 


208 Pb 82+ 


Energy / nucleon 


7 TeV 


2.759 TeV 


No. of bunches kb 


2808 


592 


Particles/bunch 


1.15X10 11 


7 x 10 7 


Transv. normalized 






omittance (la) 


3.75 /im 


1.5 [im 


RMS momentum 






spread <S 2 ) 1/2 


1.13 x 10~ 4 


1.10 x 10 -4 


Stored energy per beam 


362 MJ 


3.81 MJ 


Design luminosity C 


10 34 cm _2 S -1 


10 2Y cm- 2 s" 1 


horizontal and vertical j3* 


0.55 m 


0.5 m 



Reviews of this field can be found in Refs. [7J [H] . In 
comparison with the cross-section for inelastic hadronic 
interactions, cr/j w 8 b, the rates of ultraperipheral inter- 
actions are enormous: for example, the cross-section for 
e + e _ pair production given by the Racah formula [3] is 
of order 2 x 10 5 b. From the point of view of collider 
operation, the most important electromagnetic processes 
are the sub-classes of these interactions which remove 
ions from the beam, namely bound-free pair production 
(BFPP) and electromagnetic dissociation (EMD) [10 . 

In BFPP, the virtual photons surrounding relativistic 
ions collide and produce an e + e _ pair where, in con- 
trast to the much more frequent free pair production, the 
electron is created in an atomic shell of one of the ions. 
Schematically, the reaction between two bare nuclei with 
atomic numbers Z\, Z 2 can be written as 

Zi + Z 2 — > (Z 1 + e-) Ul/2 ,... + Z 2 + e+ (1) 

In EMD one nucleus absorbs a photon and undergoes 



2 



a transition into an excited state, typically by excitation 
of the giant dipole resonance. When this decays, it emits 
one or several nucleons. Because of the nuclear Coulomb 
barrier, the most common processes are the emission of 
one or two neutrons. For the case of 208 p D 82+ ions, the 1- 
neutron reaction (called EMD1 hereafter) can be written 
as 

208pj ;J 82+ + 208 p^S2+ > 208 + 207 p^»2+ + R Aj) 

Both BFPP and EMD change the magnetic rigidity of 
at least one of the colliding nuclei, (as usual, magnetic 
rigidity is defined as a particle's momentum p per unit 
charge Ze, or p/Ze = (Bp), where p is the bending radius 
in a magnetic field B.) If the charge of an ion changes by 
eAZ, and the number of nucleons by AA, the resulting 
rigidity can be written as (Bp) (1+6), where the fractional 
deviation 6 from the main beam with atomic number Zq 
and mass number Aq is given to a very good approxima- 
tion (neglecting increments of the mass excess) by 



X(m) 



Z (Aq + A A) 
A Q (Z + AZ) 



(1 + 6 P ) 



(3) 



Here 6 P = Ap/po is the fractional momentum deviation 
per nucleon with respect to an ion circulating on the 
central orbit of the storage ring. 

During the interactions, the transverse momentum re- 
coil is very small, so the modified ions emerge at a small 
angle to the main beam. In EMD, the recoil changes 6 P 
(see Sec. VII), while this is negligible for BFPP. However, 
ions from both processes follow dispersive orbits accord- 
ing to their magnetic rigidity and may be lost at the first 
point in the ring where the horizontal aperture A x and 
dispersion d generated locally since the IP satisfy 



d6 > A,, 



(4) 



The beam losses caused by these processes contribute 
significantly to the decay of intensity and luminosity in 
an ion collider (TUJ E] . This has been evaluated for the 
LHC [TJ Q2], where ions might collide at three IPs: the 
ATLAS experiment at IP1, ALICE at IP2 and CMS at IP5. 

The general layout of the LHC is shown in Fig. [TJ It 
consists of two concentric 27 km rings with counter ro- 
tating beams. Each ring is made of eight long straight 
sections matched to eight regular arcs through dispersion 
suppressor regions. The two rings overlap in four of the 
eight long straight sections housing colliding-beam exper- 
iments. The optics design continues to evolve through 
various versions, usually corresponding to a given lay- 
out of the collider elements, with the latest one called 
V6.503. The LHC uses superconducting magnets with 
the two apertures within the same cryostat, where some 
are cooled down to 1.9 K by liquid helium. Further de- 
tails of the LHC optics and layout can be found in Ref. |T]. 

The rate of removal of particles from the beam at the 
IPs is directly proportional to the interaction cross sec- 
tion, which in the case of BFPP takes the approximate 
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FIG. 1: The schematic layout of the LHC (the separation of 
the two rings is exaggerated). 



form [13] 



(A log 



7 cm + Bi), 



(5) 



where the electron is captured by nucleus 1 and the sum 
is taken over atomic shells i. Here 7 cm is the relativis- 
tic factor of the ions in the center-of-mass frame, and 
Ai,Bi depend only weakly on Z. For 2.76 TeV/nucleon 
208 Pb 82+ collisions in the LHC, the total cross section for 
electron capture to one of the nuclei is |13j 



281b. 



(6) 



Using Tab. [T] and Eq. ((6), the BFPP event rate is 
£<7 BFPP « 281kHz (C is the collider luminosity). 

Predictions of the cross-sections for BFPP have varied 
substantially over the years but have converged on values 
close to those found using the plane-wave Born approx- 
imation of Ref. [T3]; as these authors point out, it was 
shown earlier |14j , that higher order Coulomb corrections 
may reduce the values by a few percent . We estimate the 
uncertainty in the BFPP cross section to about 20% |13j . 

The corresponding cross section for EMD1 was calcu- 
lated with the Monte Carlo program fluka [T5l [TBI 117] . 
by simulating a large number of the events through direct 
calls to the event generator: 



96 b 



(7) 



This relies on a recent improved implementation of EMD 
effects which has been benchmarked against the RELDIS 
code [TS] (e.g., in this case reldis gives 104 b [JJ). fluka 
was also used to estimate the cross section for 2-neutron 
EMD (called EMD 2) to 



29 b, 



(8) 
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while the total EMD cross section, including all decay 
channels, is 



monitors (BLMs) to survey the losses. Finally, in Sec. IX 
we discuss possible methods of alleviation. 



226b. 



(9) 



Thus, the total cross section for particle loss by elec- 
tromagnetic processes is around 507 b, compared to the 
8 b for nuclear inelastic interactions, and is the dominant 
limit on usable luminosity lifetime for the experiments. 

Furthermore, if Eq. Q is satisfied somewhere, then the 
lost particles hit the vacuum chamber in a well-defined 
location which may lie in a superconducting magnetic el- 
ement pUCEl |20]. With the ener gy E P er ion in Tab.^J 
the total power in the BFPP beam is £a BFPP E sa 25.8 W. 
This induced heating may raise the temperature enough 
to bring the superconducting cable irreversibly over the 
critical surface in its phase diagram, the space spanned 
by magnetic field, current, and temperature. The result- 
ing departure from the superconducting state is called a 
quench. The induced Joule heating also quenches neigh- 
boring volumes so that the quench propagates. In case of 
a quench in the LHC, the beam will be dumped within 
one machine turn and the quench protection system will 
fire heaters to quickly quench a number of magnets and 
dissipate the stored magnetic energy over a larger vol- 
ume, avoiding damage to the magnets. Nevertheless, 
quenches have to be avoided by all means during collider 
operation since recovery involves the lengthy process of 
cooling the magnets down again, followed by refilling, 
ramping and "squeezing" of the beams. Downtime of the 
LHC is costly. 

BFPP has been measured in fixed target experi- 
ments [2"T1 l22l [23] - The first measurement in a col- 
lider [24] occurred during 100 GeV/nucleon 63 Cu 29+ op- 
eration of the Relativistic Heavy Ion Collider (RHIC) 
at Brookhaven National Laboratory, where the detected 
induced showers agreed within a factor 2 with FLUKA 
simulations and the location of the losses was predicted 
within 2 m. Although neither BFPP nor EMD pose any 
risk of quenches at RHIC, mutual EMD between ions of 
the colliding beams has also been measured at RHIC [25] 
and is routinely used to monitor the luminosity [25] [57] . 

The risk of inducing quenches means that it is vital 
to study, quantify and possibly alleviate the impact of 
BFPP and EMD in the LHC, in order to ensure safe op- 
eration uninterrupted by lengthy quench-recovery proce- 
dures. 

To study the loss processes, we use a three-step simu- 
lation method, consisting of optical tracking of the ions 
with modified rigidity, a Monte Carlo simulation of the 
hadronic and electromagnetic shower created as the ions 
fragment following their initial impact on the vacuum 
envelope, and a thermal network simulation of the heat 
flow inside the surrounding cryo-magnet. To illustrate 
the method we apply it to BFPP in the ALICE experi- 
ment at IP2 of the LHC, with brief comparisons with the 
other IPs in Sec. HIllVH In Sec. IVIII we treat losses from 



EMD and in Sec. VIII we describe a system of beam loss 



II. OPTICAL TRACKING 

The condition for loss of nuclei modified by BFPP or 
EMD in a localized spot, given by Eq. Q, depends thus 
on the beam optics, the aperture and the ion species. 
For 208p b 82+ iong in the LHC ^ the s caused by BFPP is 

(5 BF pp = 0.012 (using AZ = -1 in Eq. For EMD1 

(AA = -1) and EMD 2 (A A = -2) we include also the 
average momentum deviation per nucleon caused by the 
recoil, which we estimated with fluka to (8 P ) « —7.4 x 
10~ 5 for EMD1 and (5 P ) w -1.25 x 10~ 4 for EMD 2. This 
results in (5emdi = —0.0049 and <5emd2 = —0.0097 from 
Eq. @. 

Since the momentum acceptance of the LHC arcs is 
\8\ < 0.006 [1 , BFPP and EMD 2 will cause localized 
losses in the LHC while EMD1 will not. 

To illustrate the losses from BFPP and EMD, we used 
the program mad-X |28j to compute the off-momentum 
optical functions. The computation was done for both 
beams using V6.503 of the LHC optics. Fig. [2]shows the 
resulting beam envelopes for Beam 1 (circulating clock- 
wise when viewed from the top) for the nominal beam 
and the secondary beams of ions created by BFPP or 
EMD at each of the experimental IPs. 

Even if ion species other than 208 pb 82+ were used, par- 
ticles affected by BFPP would still be lost, since Eq. (§ 
and \8\ < 0.006 requires Zq > 168. The power drops 
rapidly with Z 7 according to Eq. In comparison, 

the corresponding condition at RHIC is Z > 58 so that 
the measurements in Ref. [M] were only possible with 
63 Cu 29+ beams and not the more usual 197 Au 79+ . 

Ions created by EMD1 are lost if Aq < 166, meaning 
that this might be a source of localized losses during fu- 
ture collisions between lighter ions. However, for lighter 
ions, the cross section is significantly lower |18j . 

Since BFPP is the most dangerous process in the LHC 
we focus on it here and treat EMD in Sec. IVIII 

A 208 pb 82+ ion, in the center of the bunch and follow- 
ing the ideal trajectory through the center of all mag- 
nets, enters the IP with 8 = 0. If it acquires an extra 
electron through BFPP, it exits with 8 = <5bfpp- Down- 
stream of the IP it follows a dispersive trajectory xb, 
which together with the linear optics is calculated exactly 
with MAD-X rather than using the linear approximation 
xb ~ d8. 

To write the orbit of any other 208 pb 81+ ion, we use 
subscript i to denote that the function A should be eval- 
uated at s = Si and express derivatives with respect to 
s as dA/ds — A' . Furthermore, we use a tilde to repre- 
sent chromatic optical functions for the BFPP particles 
(e.g. at s — Si the usual envelope function is Pi for the 
nominal 208 Pb 82+ beam and ^ foj . the 208 pb si+ beam )_ 

Unless stated otherwise, all optical functions refer to the 
horizontal plane. 
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FIG. 2: The simulated horizontal 6a envelopes of the nominal 208 pb 82+ beam (blue), the 208 Pb 81+ BFPP beam (green), the 
207 pb 82+ E MD1 beam (red) and the 2 06p b 82 + EMD2 beam (yellow) coming out of IP1 (top), IP2 (middle) and IP5 (bottom) 
in the LHC, with the machine aperture superimposed. The envelopes are plotted in a lighter color after the impact (indicated 
by a red star) of the central orbit on the aperture. The horizontal scales have s — at each IP in turn and the closest horizontal 
collimators are represented as black boxes at 6a to the right. The length of the collimators is not to scale for the sake of 
readability. 



With s = so at the IP, we write the position X\ and 



angle x' x of a 208 Pb 8 
s\ as 
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ion at some downstream position 



(10) 



using a standard representation of the transfer matrix in 
terms of the optical functions. The trajectory of a BFPP 
ion is thus the superposition of the central dispersive tra- 
jectory xb caused by the electron capture, a betatron 
oscillation around it, and a dispersive contribution from 
the momentum deviation within the bunch. 
The functions C\ and Si are given by 



Si 




~- [cos(/ii - Mo) + &o sin(/2i - p, ) 
Po 



/?i/3 sin (fix - /2 ), 



(11) 



—j3'/2. They are evaluated with initial conditions taken 
from the periodic on-momentum solution at the IP (/3o = 
/?o and /io = A*o)- For the off-momentum dispersion func- 
tion d, the initial condition is do = 0. 

Since the betatron amplitudes and S p are small, and 
there arc essentially no nonlinear elements in this part of 
the ring, the linear approach in Eq. ( 10 ) is accurate as 



long as the central dispersive trajectory xb is calculated 
without approximation. 

The 6 P of the BFPP particles has a Gaussian distribu- 
tion with a standard deviation given in Tab.|T] The initial 
xq and also normally distributed but we must take 

into account that the distribution of these collision points 
is not that of the incoming bunches. Assuming zero pe- 
riodic dispersion and head-on collision between identical 
Gaussian bunches at an IP, the horizontal phase space 
distribution of each bunch is 



f/3{Xo,X ) = 7^T2 ex P 



2ttct§ 



Xq + (apXQ + /3 X' ) 2 

2al 



(12) 



where fi is the off-momentum betatron phase and a = 



where Nb is the number of particles per bunch, uq = 
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\f(3oe and e is the horizontal geometric emittance. The 
luminosity density A in phase space, which corresponds 
to the distribution of the BFPP ions, is then calculated 
by integrating over the two colliding bunches, where we 
denote the coordinates in the opposing bunch at the IP 
by u and impose uq — xq as a condition for a collision to 
take place: 



X(x ,x' ) 



J fp{xp, x' o )f (x o , u' ) du' 
J fp{x , x' )fp(x , u' ) dx' Q du' dx 

go 



(13) 



This is again a Gaussian distribution, but with a smaller 
standard deviation 



Xq X(xq, x' ) dx' a dxQ 



1/2 



y/2' 



(14) 



Similarly, the standard deviation of the angular distribu- 
tion of the BFPP ions is 



e 2 



°P,0 



ft 



(15) 



which reduces to the familiar expression e//?o if ao = 
(as holds at all IPs in the LHC). We see that the distri- 
bution of collision points in space is narrower than the 
bunch distribution by a factor y/2, while the angular dis- 
tribution is similar to that of the initial bunch. The ver- 
tical coordinates y and y' are treated analogously. 

The size of the BFPP beam changes as it propagates 
through the accelerator lattice. Using Eq. ( 10 ) to propa- 



gate the beam distribution in Eq. (13 1 the standard de 



viation of the BFPP particles at s\ can be shown to be 



0X1 



sin 2 (/ii 



(16) 



This can be compared to the main beam, with 

crj = J fi\e + D\ (s){8p), where D\ is the periodic on- 

momentum dispersion. Differences come both from the 
fact that the BFPP beam envelope is oscillating with 
//, rather like a mismatched beam at injection, and the 
chromatic variation of the optical functions. 

We can deduce the impact point of the BFPP losses 
directly from Fig. [2] Once this is known, we implement 
a fast tracking of single particles from the IP to the be- 
ginning of the element where losses occur using Eq. (10 1. 



Inside this element an analytical algorithm is used to find 
the impact coordinates and momenta. 

At IP2 in the LHC, the BFPP beam converted from 
the clockwise-circulating Beam 1 is lost near the end of 
a superconducting dipole magnet in the dispersion sup- 
pressor (known as "MB" type). Tracking gives normally 
distributed losses along s, shown in Figs. [3] and [4] with a 
mean of s = 377.35 m from IP2 and a standard deviation 
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FIG. 3: The longitudinal distribution of the impact positions 
of lost 208 Pb 81 + ions from BFPP after IP2 in the LHC, as 
simulated by tracking 10 particles. The red line indicates 
the end of a superconducting dipole magnet. 
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FIG. 4: The distribution of the impact positions of lost 
20s pb 8i+ ions ( green dots ) from BFPP after IP2 in the LHC, 
as simulated by tracking 10 4 particles, shown together with 
the beam screen. 



of 65 cm. In the case of IP1 or IP5 the loss positions are 
around 418 m downstream of each IP, also in supercon- 
ducting dipole magnets, with projected rms spot sizes of 
58 cm and 74 cm respectively. The situation around the 
three IPs is therefore comparable, and we focus on IP2 
in the remainder of this paper, except where otherwise 
indicated. 

In Fig. [5] we show the on- and off-momentum (3- 
functions downstream of IP2. Around the impact point 
we have j3\ ^> $\ meaning that the BFPP beam is much 
smaller than the nominal beam. This effect concentrates 
the heat load in a smaller volume of the magnet. 



III. SHOWER SIMULATION 

The next step to determine the risk of quenches from 
the BFPP beam losses is to estimate the energy deposi- 
tion they give rise to. Defining the origin of cylindrical 
coordinates (r, 4>, z) at the center of the beam pipe at a 
magnet entrance, the heating power density P is 



P(r, 4>,z)=£ cr BFPP W(r, 4>, z), 



(17) 



where W is the average, over many possible showers, of 
the energy deposition per unit volume in the supercon- 
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FIG. 5: The on-momentum horizontal /3-function downstream 
of IP2 compared with the off-momentum 0, calculated as de- 
scribed in the text (up), and the periodic dispersion func- 
tion D shown together with the locally generated dispersion 
d (down). The dashed vertical line indicates the central loss 
position of the BFPP particles, close to the first maximum of 
the dispersion. 

ductor or other material per lost 208 pb 81+ ion. 

The quantity W, which has the unit J /(cm 3 ion), de- 
pends on the distribution in space and momentum of the 
ions lost on the beam pipe, and the cross sections for 
the many interactions that occur as the shower develops 
in the material structure of the magnet. We estimate 
W through simulations with fluka, where a 3D model 
of a main dipole has been implemented |29] as shown in 
Fig.® 

The lost particles first hit the beam screen, which is a 
racetrack-shaped chamber intended to protect the mag- 
net from synchrotron radiation. Outside is the circular 
cold bore (labelled "beam pipe" in Fig.|6|, which is sur- 
rounded by the superconducting coil, consisting of an 
inner and an outer winding. The collar is a part of the 
coil support structure and is inserted in the iron yoke. 
All these parts are heated by the hadronic and electro- 
magnetic showers. 

The fluka model of the magnet contains some sim- 
plifications: 

• The matrix of superconducting NbTi filaments and 
copper inside the coil was modeled as a homoge- 
neous body with weight fractions of the materials 
corresponding to the real coil. 

• The curvature of the magnet was neglected (al- 
though it was included in the tracking up to the im- 
pact on the beam screen). The cold mass is 15.2 m 
long and has a sagitta of 9.18 mm pQ. Simplified 
FLUKA simulations, representing the magnet as a 
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FIG. 6: The transverse cross section of the fluka model with 
the magnetic field superimposed (up), a technical drawing 1 
(down) of an LHC main dipole (middle) and the magnetic field 
in the superconducting cables (down). The rectangular boxes 
on the outside of the magnet in the fluka model represent 
beam-loss monitors as discussed in Sec. IVIIII 
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solid block (either rectangular or with a 9.18 mm 
sagitta), shows that the introduced error is around 
3-4%, which is negligible compared to the overall 
simulation error. 



1400 - 



• The magnetic field was included inside the yoke but 
omitted outside, where it cannot affect the results 
of the magnet quench evaluation. 

• Parts of the magnet far away from the coil are not 
necessary to include for the purpose of estimating 
the energy deposition in the coils, however, they 
are needed for later simulations of the beam loss 



monitors (BLMs) described in Sec. VIII 



The coordinates and momenta of the particles imping- 
ing on the inside of the beam screen, generated by the 
tracking described in Sec. |TTJ were fed as initial condi- 
tions to the fluka simulation. Some 10 4 particles were 
sufficient to keep the statistical error below 2% when the 
energy deposition in each cell of several spatial cylindrical 
grids was recorded. 

Fig. [7] shows the power density in the inner layer of 
the superconducting coil and Fig. [8] presents the energy 
deposition along the length of the cable at the hottest 
azimuth in the radial bin in the coil closest to the impact 
for different meshes. In both figures, energy deposition 
is converted to power density at the design luminosity of 
Table [T] with Eq. ([17). The maximum power density in 
the coil is P « 15.5mW/cm 3 using LHC optics V6.503, 
and does not change if the cell sizes rAr Az A<f> are re- 
duced. 

We estimate that the fluka simulation has a system- 
atic error margin of a factor 2-3, taking into account the 
uncertainties from the surface effects in the grazing in- 
cidence, the transverse momentum distribution, and the 
showering and cross sections at LHC energy. This error 
is consistent with Refs. I2~il 1501. 



IV. THERMAL NETWORK SIMULATION 
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FIG. 7: The heating power from beam losses caused by BFPP 
in the inner layer of the coil of an LHC main dipole as simu- 
lated with fluka. The power density was averaged over the 
width of the cable and is shown as a function of azimuthal 
angle (j> and longitudinal coordinate z, with z = in the be- 
ginning of the magnet. The beam loss is centered around 
z = 1206 cm and 6 ~ —3.11 rad. 
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The heat load distribution in the coil can be used to 
estimate the resulting temperature profile. Since the flux 
of lost 208 Pb 8 + ions changes only on the scale of tens of 
minutes, a steady state situation is considered, with a 
continuous deposition and evacuation of heat from the 
coil. We discuss fluctuations of the heat load in the end 
of this section. 

To better understand the heat flow, we describe briefly 
the coil geometry. The Rutherford-type cables in the coil 
(see Fig. [9]) consist of strands (28 and 36 in the inner and 
outer windings respectively), made of NbTi filaments, 
which are embedded in a copper matrix. Helium inside 
the cables occupies the space between the strands [31] but 
serves mainly to increase the heat capacity for protection 
against transient losses (32] . The cables are wrapped with 
three layers of polyimide electric insulation and a 0.5 mm- 
thick ground electric insulation is placed between the coil 



1100 1200 1300 1400 

FIG. 8: The linear power density along z caused by BFPP 
in the inner layer of the coil of the dipole magnet for varying 
radial and azimuthal binnings, where Ri <j>j stands for i radial 
bins and j azimuthal bins. In all cases, we used a longitudinal 
cell size of Az — 5 cm. The hottest bin was selected for each 
mesh. 



and the collar, which has been identified from measure- 
ments as a heat reservoir [33]. The collar is built from 
3mm-thick austenitic steel plates with a 0.1mm gap be- 
tween them filled with helium in direct contact with the 
heat exchanger through the helium in the iron yoke. The 
heat flow in steady state, schematically shown in Fig. [TO] 
is mainly limited by the heat conduction of the electric 
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insulation of the cables and the size of the helium chan- 
nels in the magnet. 

There arc several estimates of the quench limit, that is 
the smallest power density that could cause a quench. 
In Refs. [TJ |M] it is estimated as 5mW/cm 3 and 
4.5mW/cm 3 . The maximum simulated BFPP power is 
more than a factor 3 higher. These estimates are based 
on experiments |35[ 136] where a homogeneous heat depo- 
sition in a cable sample is assumed. In the case of BFPP 
however, the energy deposition is non-uniform and dis- 
tributed over the entire MB coil. 

Refs. [37J EE] describe a more detailed method of heat 
flow modeling in the coil, based on a thermal equivalent 
of an electrical network. This network model can be used 
to simulate the steady state heat flow caused by a given 
input heat source, e.g. a beam loss. It can then be in- 
ferred from the temperature map whether the magnet 
quenches or not. Since the network model simulations 
in Refs. [37, 38J show that the quench limit is highly de- 
pendent on the heat load distribution and the coil struc- 
ture, a new network model simulation, directly using the 
heat deposition map simulated with fluka, should be 
the most accurate assessment of the risk of quenches due 
to the BFPP beam losses. 

The network model takes as input a mesh of thermal 
elements, which was created from technical drawings of 
the magnet. It describes accurately the radial structure 
of the magnet from cold bore to collar, i.e., the cold bore, 
cold bore insulation, the helium channel around the cold 
bore and the coil, with the strands as the smallest ther- 
mal unit. It implements the thermal paths of the heat 
flowing from the cables through the insulation to the col- 
lar (Fig. 10). Furthermore, superfmid, normal fluid and 



polyimide tapes 50% overlapped 



gaseous helium phases are taken into account, as well as 
both heat conduction and convection, and nucleate boil- 
ing of normal fluid helium [38 . 

Other required simulation input includes the heat con- 
ductivity of the materials in the coil (calculated with 
the Cryodata software [39 ), the temperature margin 
distribution in the coil, computed with the roxie pro- 
gram |40j . and the power deposition in each strand. 
Since they are not arranged in a regular polar mesh, 
the binning used in the fluka simulation cannot be 
made to correspond exactly to the strand positions. The 
power in the strands is therefore obtained through a 
two-dimensional linear interpolation of the fluka output 
(shown in Fig. 11 1 using the cell dimensions Az = 5 cm, 
A<p = 27r/132rad, and Ar w 1.1mm (corresponding to 
14 radial bins). 

In the course of the network model simulation, the in- 
put distribution in the magnet is varied by scaling each 
value of the heat load map in Fig. [TT] up or down by 
a global scaling factor in order to determine the quench 
limit. The result of this simulation is presented in Fig. [12] 
which shows the maximum allowed input power load in 
a single strand before a quench occurs as a function of 
the magnet current for three different power load distri- 
butions caused by BFPP losses in different LHC optics 




FIG. 9: Close-up of a single cable where the current carrying 
strands, and the filaments composed of NbTi, are visible. 



Heat 




Beam 



FIG. 10: Cross section of a part of an MB magnet where the 
red arrow schematically indicates the heat flow from the coil 
to the heat exchanger. 



versions. (This is somewhat artificial since changing the 
magnet current would normally amount to a change in 
beam energy and loss distribution but the abstraction 
illuminates the physics of the magnet, considered as a 
"target".) Fig. 12 shows possible working points of the 



magnet for a given heat source distribution, and the mar- 
gin to quench at nominal beam energy can be read out 
directly. Here we discuss LHC optics V6.503 and leave 
the others for later sections. The expected heat loads are 
indicated by a horizontal lines. 



The function shown in Fig. 12 drops steeply for large 
currents, when the current itself is the limiting factor 
for quenching. In this regime, the energy is efficiently 
transported away by the superfluid helium. This is not 
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FIG. 11: The power deposition in each strand of the supercon- 
ducting coil in the 5 cm-longitudinal cell where it assumes its 
maximum value, as interpolated from the result of the fluka 
simulation. 

the case for smaller currents since, when higher temper- 
atures are tolerated, the helium loses its superfluidity. 
When not all of the energy can be transported away, the 
heating becomes the limiting factor for quenching and 
the function is approximately linear. 

As can be seen from Fig. [12] quenches are likely to 
occur at nominal operation and the luminosity for case 
V6.503 has to be decreased by 30% to go below the 
quench limit. Since the loss conditions are very simi- 
lar at IP1 and IP5, quenches can also be expected there 
at comparable levels. 

When averaging the power from BFPP losses over the 
width of the cable, the simulated quench limits for these 
loss distributions are 5-6 mW/cm 3 depending on op- 
tics version. This is in good numerical agreement with 
Refs. [UG1]) although we have assumed a higher temper- 
ature margin and lower magnetic field. Therefore, for the 
same temperature margin, our modeling of the magnet 
is more pessimistic for the BFPP heat distribution. 

All presented simulations assume steady state, where 
the heat load is constant in time. However, the real dis- 
tribution of the losses in time shows statistical fluctua- 
tions. The number of BFPP particles created in a single 
bunch crossing can be approximated by a Poisson dis- 
tribution. At nominal luminosity the expectation value 
is £cr B FPp/(fcb/rcv) = 0.043 events per crossing, where 
C and kb are given in Tab. |TJ ct B fpp by Eq. ([6]), and 
/rev = H.l kHz is the revolution frequency. The sum 
of BFPP particles created in several crossings has also 
a Poisson distribution, which can be approximated by a 
normal distribution for a large number of particles. 

To estimate if statistical fluctuations could cause 
quenches, we assume operation at 95% of the quench 
level (which might be overly optimistic) and calculate the 
probability of an increase in the number of BFPP parti- 



cles large enough to induce a quench. The network model 
simulations show that the temperature margin at this av- 
erage heat load is 0.3 K. We consider three timescales as 
in Ref. [35]: 

• t < 10 /is: During this time, corresponding to 
65 bunch crossings, 1.9 BFPP events are expected 
on average, and the heat deposited in the supercon- 
ducting cable is not yet transferred to the helium 
inside the cable. The capability of withstanding 
a quench is thus given by the cable enthalpy re- 
serve corresponding to a 0.3 K increase, which is 
calculated to 0.48 mJ/cm 3 . If the BFPP particles 
are pessimistically assumed to be lost at the same 
s- value, this corresponds to an increase in the num- 
ber of BFPP events during a 10 /is interval by more 
than a factor 1000. The probability of this occur- 
ring during 1 month of 208 pb 82+ operation at full 
luminosity (which is overly optimistic) is practically 
zero. 

• 1Q fis < t < 0.1s: The heat is transferred to the 
helium inside the cable but not out through the in- 
sulation. The enthalpy reserve is 11.7 mJ/cm 3 , cor- 
responding to an increase in the number of BFPP 
events by a factor 8 during 0.1s. The probability of 
this occurring during 1 month of 208 Pb 82+ operation 
is again vanishingly small. 

• t > 0.1s: Heat is transferred out through the ca- 
ble insulation to the heat exchanger and we con- 
sider the steady state quench limit given by the 
network model. During 0.1 s there is on average 
19068 BFPP events, while 20071 are required for a 
quench. The expected number of such fluctuations 
during one month is 4.8 x 10 -6 . 

The very small probability of a statistical fluctuation 
causing a quench justifies our steady state approach. 

A successful benchmark of the network model [3H] used 
a special heating device inserted into a cold magnet to 
produce a known heat load. Measurement and simulation 
agreed within 30% for the current where a quench occurs 
in the relevant range of heating the power. We take this a 
guide to the uncertainty of the network model simulation. 

V. SENSITIVITY ANALYSIS 

In reality, the impact points of BFPP ions vary as a 
function of imperfections such as magnet misalignments, 
field errors, and the imperfect aperture. Taking into 
account the real geometric aperture interpolated from 
measurements [41] [42] , the longitudinal loss distribution 
changes slightly, but fluka simulations show that the 
increase of peak heat load in the coils is negligible. 

We also studied the variations of the loss pattern 
caused by optical imperfections. Imperfect corrections 
resulting in a residual closed orbit caused by measured 
magnet misalignments [HJ H5J were simulated and BFPP 
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FIG. 12: Maximum allowed heat load in a single strand 
versus current of the MB dipole magnet for three different 
spatial heat load distributions caused by BFPP in different 
versions of the LHC optics: V6.500, V6.503 (current) and 
V6.503 ob (with an orbit bump as discussed in Sec. IX I. The 



horizontal lines show the heat loads expected at design lu- 
minosity as calculated by FLUKA. The required current for 
2.76 TeV/nucleon 208 Pb 82+ operation is indicated by a vertical 
black arrow. At the arrow, the tolerated power is smaller than 
the expected for V6.500 and V6.503, meaning that quenches 
can be expected during nominal operation. 



ions tracked in the resulting optics. The average longi- 
tudinal position of the impact point near IP2 varied by 
up to 2 m. A negative horizontal closed orbit, x c , in the 
BFPP impact zone causes BFPP ions to be lost further 
downstream and losses may appear in the corrector mag- 
net and beam position monitor (BPM) attached to the 
next main quadrupole. Such losses could make the BPM 
unusable. Unrealistic orbits x c < —2 mm at the impact 
point cause some BFPP ions to miss the aperture com- 
pletely at the first maximum of the dispersion function 
and continue instead to the second (see Fig. [2), where 
they hit another main dipole magnet. This is discussed 
further in Sec. IIXI 

To illustrate the sensitivity of the whole simulation 
chain with respect to optical distortions, we compared 
with V6.500, an earlier version of the LHC optics, with a 
7% smaller (3 at the impact point after IP2. However, in 
the off-momentum optics, j3 is a factor 2 smaller, reducing 
the longitudinal spot size to er\.i ~ 49 cm. The result- 
ing maximum allowed power as a function of current is 



shown in Fig. 12 This is very similar to V6.503, but the 
expected heating during nominal operation is about 25% 
higher due to the smaller spot size. 

The error sources in all simulation steps give rise to 
a large uncertainty of approximately a factor 3.9 on the 
final quench limit in the worst case. The real error is 
however expected to be much smaller than that. Our 
practical conclusion is that the heat load from BFPP 
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FIG. 13: The expected heating power from losses caused by 
BFPP at IP2 normalized by the heating power that causes a 
quench (top) and expected quench behavior at different lu- 
minosities (bottom) for various energies (labelled with the 
proton-equivalent energy), versions and configurations of the 
LHC optics. The colors are a semi-quantitative indication 
of how dangerous the losses are expected to be: Red bars 
mean that the simulated heat load is above the quench limit 
(note that operation may still be possible due to simulation 
errors), yellow that the heat load is below the quench limit 
but quenches cannot be excluded due to simulation uncertain- 
ties, while green bars can be considered as safe. The height 
of the bars indicate design luminosity. The assumed value in 
the 5 TeV configuration might be too optimistic for reasons 
of aperture, in particular for the final focusing triplets. 



beam losses are likely to be above the quench limit and 
therefore limit the luminosity of Pb-Pb collisions in the 
LHC. 



VI. COMPARISON OF OPERATING 
CONFIGURATIONS 

As the commissioning and operation of the LHC pro- 
gresses, Pb-Pb collisions will occur in a variety of config- 
urations [TJ [5] with varying beam energy, intensity, and 
other parameters such as the optical function f3* at the 
interaction point (note that 0* = 0q denned in Eq. Q). 
First ion collisions will likely take place at an energy of 
1.97 TeV/nucleon (the same magnetic rigidity as a 5 TeV 
proton beam) rather than the nominal energy of Table [T] 
where magnetic fields and excitation currents are lower, 
meaning that more heat can be absorbed before a quench 
takes place (see Fig. 12 1. On top of that, both the lumi- 



nosity and the energy deposited per lost BFPP ion will 
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be lower. 

We repeated the tracking and fluka simulation for 
this case, assuming /3q — 0.5 m at IP2 (which is over- 
optimistic for reasons of aperture but a useful compar- 
ison). The resulting distribution of the heat deposition 
is very similar to the nominal case apart from a global 
scaling factor of 0.63. Thus there is no need to redo 
the thermal network simulation — the maximum tolerated 
heat load can be read out directly from the curve 6.503 



in Fig. 12 at the appropriate current of 8.46 kA. The er- 
ror introduced by omitting the last simulation step is less 
than 10% if we use as a benchmark V6.500, for which the 
full simulation chain was carried out. This gives an ex- 
pected heat load with a 1.97 TeV/nucleon 208 Pb 82+ beam 
at 46% of the quench level, which means that quenches 
are unlikely to take place but still within the simulation 
uncertainty. 

We have performed similar tracking and fluka simu- 
lations for various values of po that may be used for col- 
lisions. Increasing /3 reduces the luminosity (C cx Pq 1 ) 
and the event rate for BFPP (see Eq. ([l7)). Further- 
more, in each optics version /3 and d at the impact point 
after IP2 are slightly different, thus causing variations in 
the loss distributions. The heating in the magnet scales 
therefore only approximately as P(x, y, z) oc Pq 1 ■ 

The results of the quench performance for BFPP at 
IP2 for all optics versions are summarized in Fig. |13| 
The bottom part shows expected luminosity limitations, 
where we have assumed that the spatial heat load dis- 
tribution in a certain distribution stays constant when 
the luminosity varies. This is strictly true only when 
the luminosity is changed through a decrease in bunch 
population or number of bunches. 

In Fig. [13] we also show the result for the nominal op- 
tics but with an orbit bump introduced to reduce the 
maximum energy deposition. This is discussed in detail 
in Sec. IXj 

No heavy ion collisions are planned to take place at 
injection energy. 



VII. LOSSES FROM EMD IN THE LHC 

Since the ions affected by EMD1 stay within the mo- 
mentum acceptance of the LHC ring, they are intercepted 
by the momentum collimation system (see Fig. [2]) . 

For EMD 2 ions, we consider 8 p as the sum of two ran- 
dom variables: The natural momentum deviation in the 
bunch with standard deviation a « 1.1 X 10~ 4 [Q and the 
recoil caused by EMD, which was simulated with fluka 
(see Fig. 14 1. The resulting distribution is approximately 
Gaussian with mean /i = —1.25 x 10~ 4 and standard de- 
viation a = 2.9 x 10~ 4 . 

We used Eq. (10 1, with xb replaced by the central 



EMD2 trajectory, to track particles from IP2 with this 
^-distribution and a phase space distribution given by 



Eq. (13). The resulting loss map is shown in Fig. 15 
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FIG. 14: The distribution of 5 P caused by EMD2 as simulated 
by fluka. The distribution has (S p ) £s —1.25 x 1CT 4 and 
V(Vp ~ <<5p» 2 > « 2.6 x 10- 4 
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FIG. 15: The distribution of 206 Pb 82+ losses after IP2. The 
two bottom plots show zooms on the largest peaks on the 
upper plot. Only the loss at s ~ 3125 m risks to induce 
quenches, because of the small projected spot size. 



tions, where there is no risk of quenches due to the small 
intensity and the wide spot sizes, while the other half is 
lost at s w 3125 m towards the end of a drift, 70 cm in 
front of an orbit corrector magnet, called MCBC. Since 
the losses occur at a place where the aperture is steeply 
decreasing, they are distributed over a very small lon- 
gitudinal distance with a standard deviation of 2.9 cm. 
Therefore, the heat density per particle is correspond- 
ingly higher than in the case of BFPP. In Fig 



16 



Around half of the losses is distributed over several loca- 



we 

show a close-up of the impact of the 6-cr envelope of the 
EMD2 beam from IP2. 

The quench limit of the corrector is not well known. It 
is a function of the current and magnetic field, which are 
highly dependent on the closed orbit and the global orbit 
correction scheme. Furthermore, no network model of 
the corrector exists, so detailed simulations with FLUKA 
and the network model were not carried out. Assuming 
its quench limit to be similar to the dipole's, we conclude 
that there is a risk of quenching from the small spot size. 
Since the beam is automatically dumped if the corrector 
quenches, this can not be allowed to happen. 

However, it is possible to exclude this magnet from the 
global orbit correction. This means that it will not be 
operating during heavy ion runs, excluding quenches at 
the price of possible orbit distortions. Therefore losses 
from EMD2 are less serious than from BFPP, since a 
dipole is necessary for machine operation. During the 
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FIG. 16: The impact of the 6a envelope of the 206 Pb 82 + ions 
created by EMD2 at IP2 at an aperture restriction. The pro- 
jected longitudinal rms spot size on the beam screen is only 
2.9 cm. The loss takes place in a drift section and the be- 
ginning of next magnet, a corrector, is indicated by a dashed 
vertical line. 




FIG. 17: Interior part of an ionization chamber used as a 
beam loss monitor in the LHC. Taken from Ref. 43 . 

could well be different depending on the type of particle 
lost and the showers it gives rise to. Accordingly, there 
are two important problems related to the monitoring of 
ion beam losses from the collisions: 

1 . To determine the BLM threshold signals for Pb ion 
losses. 



first ion runs at 1.97 TeV/nucleon, the corrector is much 
less likely to quench than at top energy, for the reasons 
explained in Sec. |VI| 

Further losses from EMD2 take place 433.7 m down- 
stream of IP1 and IP5, in a drift section only 30 cm 
in front of a main quadrupolc. In this case the pro- 
jected rms beam sizes of the impinging ions are 30 cm 
for both IP1 and IP5. In spite of this small size, the 
risk of quenches is small for several reasons. The quench 
limit of a quadrupole is around 40% higher than that of 
a dipole [35] and the cross for EMD2 section is about 
10 times smaller than for BFPP. Furthermore, during 
the measurements described in Ref. [38 , it was found 
that the tolerable heat load is approximately 50% higher 
in the ends of the magnet, where the magnetic field is 
lower and the cooling more efficient. However, to quanti- 
tatively determine the margin to quench, network model 
simulations should be carried out. 



VIII. MONITORING LOSSES 

Since beam losses caused by BFPP or EMD might in- 
duce quenches, it is vital to survey these losses while 
beams are colliding and to make sure that the beam is 
extracted quickly to a dump before a quench can occur. 
The LHC's beam-loss monitor (BLM) system has been 
designed to detect losses around the ring during proton 
operation [43 , 44 . The LHC BLMs are 50 cm long ioniza- 
tion chambers filled with nitrogen that detect secondary 
charged shower particles outside the magnet cryostat. In 
order to minimize the drift time of the ions and elec- 
trons set free in the chamber, there is a stack of parallel 
aluminum plates inside the cylinder with alternating po- 
larities and 5 mm spacing. A view of the interior part of 



an ionization chamber is shown in Fig. 17 



The locations of the BLMs and the thresholds for trig- 
gering the beam dump system have been determined for 
protons through simulations [45j 06] . The BLM thresh- 
old signal should correspond to a a certain level of power 
deposition in the coil and the ratio of these two quantities 



2. To verify that the placement of BLMs provides ad- 
equate coverage of the loss patterns. 

To investigate the first problem, we simulated the de- 
velopment of the showers generated by particle losses, 
both from 208 pb 82+ ions and protons, in an LHC dipole 
magnet with fluka; the geometry was as described in 
Sec. |III| and Fig. [6] The BLMs are schematically modeled 
as thin rectangular iron boxes filled with nitrogen, placed 
outside the MB cryostat. This simplification is consistent 
with simulations done for protons |43] and can be consid- 
ered reasonable since we are primarily interested in the 
ratio between heat deposition from heavy-ion and proton 
losses rather than the absolute signal from the BLM. 

For each particle species, a generic beam loss was rep- 
resented by a "pencil" beam: a monoenergetic beam of 
particles, all hitting the vacuum chamber at the same 
point (in the horizontal plane, 1 cm from the entrance of 
the magnet) at the same angle of incidence (0.5 mrad). 
An arbitrary loss can be modeled as a superposition of 
pencil beams. Simulations were done with 208 pb 82+ ions 
and protons at 2.76 TeV/nucleon and with protons at 
7 TeV. 

During the simulation, the energy deposition was 
scored in the beam screen, in the superconducting coil 
(using the mesh described in Sec. Ill I and in the BLMs. 



Since the BLM signal is proportional to the ionization en- 
ergy loss in its gas volume and ionization is the dominant 
energy loss process for the low energy charged secondaries 
that emerge outside the cryostat, this is a fair approxi- 
mation. The resulting energy deposition profiles along 
the hottest superconducting cable in the MB coil and in 



the closest BLM are shown in Fig. 18 The ratio between 
the heat deposited in the coils and the energy deposition 
in the BLMs is similar for 208 pb 82+ ions and protons at 
the same energy per nucleon, while 7 TeV protons have 
an even lower ratio between BLM signal and heating of 
the coil than 2.76 TeV/nucleon 208 Pb 82+ ions. 

This similarity comes from the fact that the particles 
causing ionization of the gas in the BLMs are not the 
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Pb 82+ ions and protons lost directly from the beam, 
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FIG. 18: The energy deposition in the hottest wire in the coil 
and in the N2 gas inside the ionization chamber for Pb ions 
and protons. The three cases are scaled to equivalent total 
energy of incoming particles. 
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FIG. 19: The energy deposition for the 208 pb 82 + ions and 
protons at LHC energy in the beam screen, as a function 
of the longitudinal coordinate in the magnet, averaged over 
0.1 x 0.1mm 2 transversally and 10 mm longitudinally. 



but instead low-energy secondary particles, mainly elec- 
trons, created in the hadronic shower. The energy de- 
position resulting from the hadronic shower is very simi- 
lar for heavy ions and protons |47j . even though the nu- 
clear interaction lengths are very different (0.8 cm for 
2.76 TeV/nucleon 208 Pb 82+ and 25 cm for 7 TeV protons 
on a Cu target according to fluka simulations). This 
comes from the fact that when an ion traverses a tar- 
get, the nucleus splits up into smaller fragments through 
electromagnetic dissociation and nuclear interactions in 
several steps. Once totally fragmented it gives rise to 
a similar shower profile as independent nucleons. The 
energy deposition from the first few interactions is far 



exceeded by deposition by the low energy secondary par- 
ticles created later in the shower. 

If the energy per nucleon increases, so do the cross sec- 
tions, and more energy is deposited closer to the impact. 
This explains why 82 protons at 7 TeV have a somewhat 
higher energy deposition in the coil than the 208 protons 
at 2.76 TeV. 

The only location where a significant difference in en- 
ergy deposition between 208 pb 82+ ions and protons can 
be seen is in the beginning of the shower, close to the cen- 
tral core. Here the difference in ionization energy loss, 
given approximately by the well-known Z 2 dependence 
in the Bethe-Bloch formula |H], is clearly visible and the 
ions cause a much higher energy deposition. Thanks to 
the small impact angle, this effect is only visible in the 
beam screen, as shown in Fig. [19] 

Based on the similar ratio between energy deposition 
in the coil and in the BLMs outside the cryostat, we 
conclude that we can use the same beam dump thresholds 
for 208 Pb 82+ ions and protons in the LHC. This result 
applies to all mechanisms for ion beam losses in the LHC, 
not only those discussed in this paper. 

Suitable positions for BLMs were determined from the 
optical studies illustrated in Fig. [2j Studies were per- 
formed both for the nominal parameters and a configu- 
ration known as the "early ion scheme" optics (which has 
a higher (3* = 1.0 m, see Chap. 21 in Ref. [1J) and for the 
two beams circulating in opposite directions. 

Because of the uncertainty of the impact position, de- 
scribed in Sec. [V] and the fact that the loss peaks are 
narrow and localized, the LHC's machine protection sys- 
tem needs a very tight coverage with BLMs. Based on 



the energy deposition profile in the BLM gas (see Fig. 18 1 
a 1.5 m spacing between chambers has been assumed for 
both beams to ensure full detection and localization of 
losses. 

The BLMs foreseen for proton operation are mounted 
on all quadrupoles in the arcs and dispersion suppres- 
sors [35] . Extra chambers for monitoring ion losses from 
BFPP have been requested in all loss locations that are 
not already covered. Since a fraction of the losses might 
escape the first loss location and instead hit the aperture 
close to the second peak of the dispersion, this position 
will also be monitored. 

No extra BLMs are needed for losses from EMD2, 
since the loss locations are already covered by the de- 
fault scheme. 



IX. ALLEVIATION OF LOSSES 

Ideally, the optics and aperture of a heavy-ion collider 
would be such that the transformed ions stay within the 
acceptance of the ring and arc cleaned by the collimators. 
But this may conflict with other design constraints like 
the overwhelming cost advantage of keeping the aperture 
small. In the case of the LHC, quenches caused by lost 
BFPP and EMD2 ions are likely and methods to avoid 
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them have to be found. Several methods using specialized 
hardware are possible, the most promising one being the 
installation of extra collimators around each IP where the 
dispersion starts to rise in the cold part of the machine [1 , 
H9"j . This could be a very efficient solution but requires 
extra hardware to be installed, as well as displacing some 
magnets to make room in the tunnel and cannot be done 
before the first heavy-ion physics runs. 

For these runs at least, other methods to mitigate the 
effects of BFPP and EMD2 using existing hardware will 
be needed. We cannot influence the quench limit of the 
dipole magnet, nor er BFPP , which leaves only L and W 
as the free parameters (see Eq. (17l). It has been sug- 
gested that, for certain conditions on the filling time of 
the LHC, a higher integrated luminosity may be reached 
if the peak luminosity is reduced [50 by varying (3* dur- 
ing a fill as a function of the remaining intensity. How- 
ever this only pays off if the single-bunch intensity can be 
raised. Otherwise the only remaining option is to reduce 
W by varying the distribution of the impacting losses. 

This might be achieved by manipulating the 
quadrupoles between the IP and the impact point to in- 
crease (3 or d, spreading out the spot of heat deposition. 
However, strict requirements on the betatron phase ad- 
vance around the LHC [51] . matching conditions on the 
periodic (3 and dispersion functions at the ends of the dis- 
persion suppressors and constraints on magnet strengths 
limit the scope of this approach. A gain in spot size of 
at best 20% was achieved using this method. 

Another alleviation method might be to adjust the or- 
bit or optics to move the impacts to more favorable po- 
sitions. The horizontal orbit of the particles affected by 
BFPP and EMD2 emerging from IP2 in a perfect lat- 
tice, as calculated by mad-x, is shown in the top part of 
Fig. [20] together with the main circulating beam. This 
is a close-up of the first 700 m in Fig. [2] The dispersive 
trajectories oscillate with d, and the BFPP particles are 
lost very close to its first maximum, which we call s±, 
while the EMD2 beam continues further downstream. 



The green BFPP envelope in Fig. 20 is significantly 
wider at the second maximum s 2 of d, thanks to the larger 
value of (3 (see Fig. [5]) . A localized closed orbit bump in- 
troduced around s± displaces xb towards the center of 
the beam pipe and allows some of these particles to es- 
cape downstream to S2- There, the lost ions are diluted 
over a larger volume decreasing the maximum power de- 
position. 

At the same time, the envelope of the EMD2 ions 
moves towards the outside of the beam pipe at s% pro- 
voking losses there. On balance this is beneficial because 
the aperture is constant in the neighborhood of s 1; so 
the losses are spread out much more than at the later 
loss point shown in Fig. 16 This new loss position for 



the EMD2 ions is the same as for BFPP ions with a flat 
orbit but since the EMD2 ion flux is almost a factor 10 
lower than for BFPP, the heating at this position is kept 
low enough to avoid quenching the dipole (see Fig. 12 1. 
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FIG. 20: Horizontal physical aperture and 6a envelopes after 
IP2: with the usual flat orbit (top), with an orbit bump (mid- 
dle), and with decreased dispersion (bottom). The beams 
are the main circulating 208 pb 82 + beam (blue), the 208 Pb 81 + 
BFPP beam (green), and the 206 Pb 82 + EMD2 (yellow). The 
vertical dashed lines show the locations of horizontal corrector 
magnets and those indicated in red are active (in this case, 
correctors called "MCBC" and "MCB"). 



rectors attached to the main quadrupoles in the LHC. 
The orbit bump amplitude required was calculated to 
keep 95% of the BFPP beam, with the beam size given 
by Eq. (16 1, inside the aperture at s\. The resulting 
beam envelope with the bump active is shown in Fig. |20[ 



The off- momentum orbit can be displaced by orbit cor- 



The orbit of the circulating beam is displaced by up 
to 3.8 mm in this configuration but it is clear that the 
6cr envelope is still far from the aperture. Apertures in 
the LHC are conventionally characterized using a quan- 
tity rii, defined as the maximum acceptable primary col- 
limator opening, in units of beam a that still provides 
a protection of the mechanical aperture against losses 
from the secondary beam halo [1]; the complete defini- 
tion incorporates a variety of tolerances that we shall not 
enumerate here. The bump reduces ri\ for the circulat- 
ing beam, from 30 to 21. This still provides sufficient 
aperture margin. The three corrector magnets used to 
make the bump would operate at between 9% and 26% 
of their maximum strength, leaving a comfortable mar- 
gin for their function in the global orbit correction. The 
/3-beating due to the bump is around 0.3%, which is ac- 
ceptably small pp. 

We have repeated the simulation chain of tracking, 
fluka, and network model for the BFPP ions with the 
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orbit bump included. The spot size at the new im- 
pact position, in another superconducting dipole mag- 
net, is 220 cm and the maximum heat load from FLUKA 
is 4.5mW/cm 3 , about a factor 3.5 less than without the 
bump. The result of network model simulation is in- 
cluded in Fig. [12] The expected heat distribution at S2 
with the orbit bump activated (indicated by a horizontal 
line) is predicted to be a factor 2.25 below the quench 
limit, still not completely safe. 

Operationally, the bump amplitude will have to be 
fine-tuned around the theoretical value to compensate 
the local closed-orbit and misalignments of the beam 
pipe. The BL Ms pla ced around each loss location (as de- 
scribed in Sec. VIII I will be an essential tool to monitor 



where a;, is given by Eq. (10) as a function of the initial 



how the BFPP losses move from the first to the second 
maximum of the dispersion. 

The orbit bump will be set up at low luminosity, in 
order to avoid potential quenches, by means of a van der 
Meer scan |52j . in which the beam orbits are scanned ver- 
tically across each other at the IP. One could also imag- 
ine tuning the bump using a higher value of (3* at the IP. 
However, this is less reliable since the optical functions 
between the IP and the impact point will change with (3* . 
Another option is to introduce the bump at lower beam 
intensity in an earlier fill and rely on reproducibility. 

From the point of view of machine protection, we must 
consider the possibility that one of the three correctors 
used in the bump might quench, leaving the bump open 
and possibly causing damage. Two of them (of type 
MCBC) are directly connected to the beam abort sys- 
tem which would dump the beam immediately. In case 
of a quench of the third corrector (of type MCB), no au- 
tomatic beam dump occurs. The resulting global orbit 
distortion will, if it is small enough, be corrected by the 
orbit feedback system, in which case the BFPP beam 
may quench a magnet. If the distortion is too large, the 
resulting beam losses will trigger a beam abort via the 
BLMs (and possibly also the beam position monitors). 
In all cases, the quench protection system should protect 
the magnets from physical damage. 

If moving the losses to another location with a larger 
spot size was insufficient, one might spread out the losses 
with several orbit bumps, which allow fractions of the 
BFPP beam past each maximum of d. By tuning the 
orbit bumps the losses can be spread over n dispersion 
maxima by n — 1 orbit bumps using n + 1 correctors in 
such a way that, ideally, a fraction l/n is lost at each im- 
pact point. This decreases the maximum heating power 
in a single element by 1 jn and can bring it below the 
quench limit if n is made large enough. However, since 
BLMs have to be used when the bumps are tuned, the 
precision of the achieved loss distribution is limited. 

To determine the required bump amplitudes A; at each 
impact location Si we consider the initial phase space at 
the IP. The particles lost at a location with horizontal 
aperture Ai satisfy (for i £ [1, n — 1]) 



conditions at the IP. These inequalities define regions Ri 
in the initial phase space, which vary with Aj. The frac- 
tion Fi of particles lost at Si is the integral distribution 
function over the phase space area outside the aperture 
limitation ( 18 1 and not outside any previous aperture 
limitation: 

Fi = / X(x ,x' )fsp(S p ) dx dx' dS p 

Ji?;n(fi;5;u^...ufl c _ x ) 

(19) 

where R c denotes the complement of region R and f$ p is 
the assumed Gaussian distribution function of 5 p . 
The Ai can then be determined by requiring 

Fi = l/n, Vz (20) 



and solving Eqs. (19 1 recursively, starting at i — 1, which 



can be solved analytically to yield 



Ai = \/2(j\ 1 erf 



n - 2 



xb.i - A 1 . 



(21) 



A, > At 



(18) 



Here cta.i is given by Eq. ( 16 1 . Numerical integration and 
solution has to be used fornighcr values of i. 

The method of one or several orbit bumps only works 
as long as the displacement of the closed orbit is kept 
within acceptable limits and the corrector magnets do 
not operate too close to their maximum strength. This 
is not the case at IP1 and IP5, where the required bump 
amplitude for BFPP particles to avoid the first maximum 
of d is 8-9 mm. The method will not work there unless 
other changes are made to the optics. 

Instead of an orbit bump, another method to move 
losses to S2 would be to tune the quadrupoles after IP2 in 
order to decrease the dispersion at s\. mad-x was used to 
rematch IP2 with all constraints mentioned in Ref. |51j 
and beam envelopes for one solution are shown in the 
third plot of Fig. [20] The longitudinal size of the spot at 
S2 is 179 cm, meaning a higher heating power than with 
the orbit bump. So quenches cannot be excluded with 
this method either. 

The required change in phase advance in the insertion 
is 0.05 x 2-7T, which can be compensated in the insertion 
in IP4 containing the RF system J5TJ to keep the overall 
tune in the machine constant. This method has the ad- 
vantage of not displacing the circulating beam but does 
not reduce the far-downstream losses from EMD2. Given 
the numerous other constraints on the LHC optics, it is 
difficult to envisage applying this method at IP1, IP2, 
and IP5 simultaneously. 



X. CONCLUSIONS 

Electromagnetic interactions such as bound-free pair 
production and electromagnetic dissociation in ultrape- 
ripheral nuclear collisions at the LHC, modify the charge 
and mass of beam ions. These particles follow disper- 
sive orbits until they are lost in locations determined by 
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the machine optics, aperture and the magnetic rigidity 
of the ions. When sufficiently localized they can heat su- 
perconducting magnetic elements enough to make them 
quench. 

We have presented the first fully-integrated simulation 
chain of beam tracking, shower simulation and a ther- 
mal network model to evaluate the heat-flow and quench 
behavior of the superconducting coils immersed in su- 
pcrfluid liquid helium. This simulation has been applied 
to the most critical loss mechanism, BFPP, occurring in 
the three heavy-ion collision points in the LHC. Heat de- 
position caused by 208 pb 81+ ions in main dipoles down- 
stream from IP1, IP2, and IP5 is expected to be 40% 
above the quench limit, while 207 Pb 82+ from EMD1 stay 
within the acceptance of the arc and are cleaned by the 
collimation system. Furthermore, depending on its exci- 
tation level, there is some risk from quenches of a correc- 
tor magnet downstream of IP2 by 206 pb 82+ ions created 
through EMD2. 

To avoid quenches, an efficient beam- loss monitor sys- 
tem to detect these losses is needed. Shower simulations 
of the relations between the signals in the loss monitors 
and the energy deposition in superconducting coils have 
been carried out for 208 Pb 82+ and proton beams. These 
demonstrate that, thanks to nuclear fragmentation, the 
beam abort system can be set to trigger at the same sig- 
nal level for both beams. Additional monitors have been 
installed in critical locations for heavy-ion losses in the 
LHC. 



Finally, we have investigated methods to alleviate the 
losses caused by BFPP and EMD2. We can move them 
to different locations in the machine, where the losses are 
spread out over a larger distance, lowering the energy de- 
position to around 45% of the quench limit. Approaches 
using orbit correctors to create a local orbit bump or 
quadrupolc tuning to decrease the dispersion have been 
illustrated. These methods have the advantage of not 
requiring new hardware in the LHC but are not easily 
applied to all interaction points and do not provide suffi- 
cient safety margins against quenches. An efficient solu- 
tion is likely to require new hardware such as additional 
collimators or masks in the dispersion suppressor sec- 
tions. 
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